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Abstract. Quasi-periodic oscillations (QPOs) have been discovered in 
the X-ray emission of many neutron stars and black holes. The QPOs 
with frequencies ^ 300 Hz are thought to be produced near the surfaces 
of neutron stars and the event horizons of black holes. I first summarize 
some of the most important properties of the QPOs seen in neutron star 
and black hole systems. I then review some of the models that have been 
proposed and compare them with observational data. Finally, I describe 
how these QPOs can be used to determine the properties of dense matter 
and strong gravitational fields. 



1. Introduction 

It is an honor to be able to contribute to this volume memorializing Jan van 
Paradijs. Jan made important scientific contributions to many areas of astro- 
physics, but especially to our understanding of thermonuclear X-ray bursts and 
the properties of low- mass X-ray binary systems (LMXBs). It is therefore par- 
ticularly fitting to include a few words about Jan's scientific contributions in 
this review, which is concerned with X-ray oscillations in LMXBs. 

Jan was one of the first to recognize that careful study of X-ray bursts could 
provide quantitative information about the properties of neutron star matter and 
strong gravitational fields. He realized that the Rossi X-Ray Timing Explorer 
(RXTE) would make possible much more sensitive and detailed studies of these 
questions. Jan led a collaboration that anticipated the detection of periodic 
brightness oscillations during thermonuclear X-ray bursts and 10-60 Hz QPOs 
in the atoll sources. Even before RXTE was launched, he and his collaborators 
requested and were awarded observing time to search for these two phenomena. 
Both were discovered soon after RXTE was launched and are among the most 
important discoveries made with it. Jan's team was somewhat unlucky in not 
catching any bursts with detectable oscillations early in the mission, but they 
discovered kilohertz QPOs in many neutron star X-ray sources. As I describe 
below, these and other discoveries made with RXTE have opened a new era in 
the study of dense matter and strong gravitational fields. I continue to miss Jan 
ELS £1 friend as well scientific colleague. 

Neutron stars and black holes are important cosmic laboratories for study- 
ing fundamental questions in physics and astronomy, especially the properties 
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of dense matter and strong gravitational fields. The discovery using RXTE 
of oscillations in the X-ray emission of accreting neutron stars and black holes 
with frequencies comparable to the dynamical frequencies near these objects has 
provided important new tools for studying these questions. 

The neutron stars in LMXBs have proved to be particularly valuable sys- 
tems for investigating the innermost parts of accretion disks, gas dynamics and 
radiation transport in strong radiation and gravitational fields, and the prop- 
erties of dense matter, because the magnetic fields of many of these stars are 
relatively weak (~ 10 -10 10 G) but not negligible. Magnetic fields of this size are 
weak enough to allow at least a fraction of the accreting gas to remain in orbit 
as it moves into the strong gravitational and radiation fields of these stars, but 
strong enough to produce anisotropic X-ray emission, potentially allowing the 
spin rates of these stars to be determined. The "clockwork" of these systems is 
remarkably regular and systematic. Observations using RXTE are providing a 
wealth of diagnostic information about strong gravitational fields and the mag- 
netic fields, spin rates, masses, and radii of these stars (for previous reviews, see 
Lamb, Miller, & Psaltis 1998a, 1998b; van der Klis 1998, 2000; Psaltis 2000). 

QPOs have also been discovered in the X-ray emission of accreting black 
holes in LMXBs. Study of this variability is providing increasingly detailed 
information about these objects. The recent discovery that GRO J1655— 40 
and GRS 1915+105 both have two high-frequency QPOs (Remillard et al. 2001; 
Strohmayer 2001c, 2001d) is the first evidence that high-frequency QPOs also 
occur in pairs in black hole systems and the strongest evidence to date that the 
black hole GRO J1655— 40 has significant angular momentum. 

The discovery with RXTE of X-ray variability on dynamical time scales 
in black holes and neutron stars has revolutionized the study of accretion and 
X-ray emission by these objects. Before the launch of RXTE, Newtonian mod- 
els were often used; after the discovery of very high frequency oscillations with 
RXTE, strong-field general relativistic calculations have become the norm, to 
allow quantitative comparisons of models with the RXTE data. The highest- 
frequency QPOs are thought to come from regions where the gravitational field 
is very strong. Study of these oscillations may provide the best near-term op- 
portunity to observe and measure uniquely general relativistic effects, such as 
gravitomagnetic 1 (Lense-Thirring) precession and the absence of stable circu- 
lar orbits around sufficiently compact relativistic objects, and to determine the 
properties of the dense matter in neutron stars. RXTE continues to provide a 
wealth of new information on conditions close to the surfaces of neutron stars 
and the event horizons of black holes and to stimulate imaginative new thinking 
about the physical processes that dominate in these extreme conditions. 

In §2 I summarize some of the most important properties of the X-ray os- 
cillations seen in neutron star and black hole systems. In §3 I review some of the 
mechanisms that have been proposed for generating these oscillations, compare 
them with some of the data on neutron star and black hole QPOs, and discuss 
the similarities and differences of QPO generation in neutron star and black hole 



lr The term "gravitomagnetic" refers to the close analogy between the motion of a test particle 
in the spacetime near a circulating mass current and the motion of a charged particle in the 
magnetic field near a circulating electrical current. No magnetic field is involved. 
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systems. In §4 I describe how measurements of the highest-frequency QPOs can 
be used to determine the properties of dense matter and strong gravitational 
fields. In §5 I conclude with a few remarks about the progress we have made 
in understanding black hole and neutron star QPOs and using them as tools to 
better understand the properties of strong gravitational fields and dense matter. 
There is space here to review only a fraction of the large body of work on these 
problems. I apologize to those whose work is not adequately discussed. 



2. Properties of Neutron Star and Black Hole QPOs 

This section summarizes several important properties of the oscillations that 
have been detected in the X-ray emission of accreting neutron stars and black 
holes. More detailed reviews of the observations may be found in van der Klis 
(1998, 2000). 

2.1. Brightness Oscillations in Neutron Stars 

Three classes of X-ray brightness oscillations have been detected in neutron stars 
in LMXBs: nearly periodic oscillations that, with one exception, have been 
detected only during thermonuclear X-ray bursts; pairs of kilohertz QPOs in 
accretion-powered emission; and several lower-frequency QPOs, also in accretion- 
powered emission. A key feature of the QPOs and peaked noise components 
seen in neutron stars is that the frequencies of all of them (except the recently 
identified hectohertz peaked noise component) vary significantly and are posi- 
tively correlated with each other and with the inferred X-ray flux (see, e.g., van 
Straaten et al. 2001). 

X-ray burst oscillations. — X-ray flux oscillations with frequencies ranging 
from ~ 330 Hz to ~ 590 Hz have been detected during some of the thermonu- 
clear X-ray bursts of ten neutron stars (see Strohmayer 2001a, 2001b). These 
oscillations remain nearly sinusoidal and relatively stable (Ais/v ~ 0.01) even 
when their peak-to-peak amplitudes are ~ 70% and the luminosity of the star 
approaches the Eddington critical luminosity at which the outward force of the 
radiation balances the inward force of gravity. This indicates that they are not 
caused by oscillations of the stellar surface or the accretion disk, because such 
oscillations would become nonlinear at high amplitudes and would also be af- 
fected by strong radiation forces. Their stability is easier to understand if their 
frequencies are related to the spin of the star, which is very stable. 

The X-ray spectra, waveforms, and long-term stability of these oscillations 
indicate that they are generated by rotation of a brighter region, or two antipo- 
dal brighter regions, of the stellar surface (Strohmayer et al. 1996, 1997, 1998; 
Bildsten 1998; Miller & Lamb 1998; Markwardt, Strohmayer, & Swank 1999; 
Miller 1999; Strohmayer & Markwardt 1999; Strohmayer 1999, 2001a, 2001b; 
Weinberg, Miller, & Lamb 2001). If so, the frequency v osc of the strongest os- 
cillation is close to the stellar spin frequency v s or its second harmonic. The 
frequencies of the burst oscillations are comparable to the spin frequencies pre- 
dicted by the recycling hypothesis for millisecond rotation-powered pulsars and 
the theory of disk accretion by weakly magnetic neutron stars (see Ghosh & 
Lamb 1992 and references therein). 
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Several properties of the burst oscillations are not yet understood. These 
properties include the relatively large fractional changes in frequency (up to 
~ 1.3%) seen in some sources (Galloway et al. 2000), the simultaneous presence 
of several distinct frequencies (Miller 2000a; Galloway et al. 2000; Chakrabarty 
2000), phase shifts during some bursts that are ~ 10-7T (Miller 1999b), and very 
low upper limits on the amplitudes of higher harmonics of the basic oscillation 
frequency in some sources (Strohmayer 2001a, 2001b). 

A strong (amplitude ~ 50%) brightness oscillation has recently been found 
in data taken during a 1997 thermonuclear X-ray burst from SAX J1808— 369 
(in't Zand et al. 2001). This source is an accretion-powered X-ray pulsar (Wij- 
nands &: van der Klis 1998) in a low-mass binary system (Chakrabarty & Morgan 
1998). The frequency of the burst oscillation is within a few Hertz of the 401 Hz 
spin frequency inferred from the periodic modulation of its accretion-powered 
emission. The burst oscillation could be produced by a mechanism different from 
the mechanism that produces the burst oscillations seen in other neutron-star 
LMXBs, but this seems unlikely. Thus, although the detailed physics of the 
burst oscillations is not yet understood, the evidence that they are generated by 
the spin of the neutron star is strong. 

Kilohertz QPOs. — QPOs with frequencies between ~ 500 Hz and ~ 1300 Hz 
have now been detected in the accretion-powered emission of more than twenty 
neutron stars in LMXBs (see van der Klis 2000). They usually appear as a 
pair of prominent, narrow peaks in power spectra. The frequencies v\ and v-i of 
the lower and upper kilohertz QPOs move upward and downward together and 
can increase and decrease by as much as a factor of two within a few hundred 
seconds while remaining highly coherent (f/FWHM ~ 100; see van der Klis 
1995 and Wijnands et al. 1998), indicating that the QPO mechanism is highly 
tunable. Their frequencies are strongly positively correlated with the X-ray 
flux. In some sources, the frequency separation Av = v<i — v\ between the two 
kilohertz QPOs decreases systematically by 30-100 Hz, depending on the source, 
as V2 increases by a much larger amount (Sco X-l: van der Klis et al. 1997; 
4U 1608-52: Mendez et al. 1998; 4U 1735-44: Ford et al. 1998; 4U 1728-34: 
Mendez & van der Klis 1999; see also Psaltis et al. 1998b). The kilohertz QPOs 
have rms amplitudes in the 2-60 keV band of the RXTE Proportional Counter 
Array ranging from <; 1% up to ~ 15% and quality factors Q = I//FWHM as 
large as ~ 200. Sidebands of the lower kilohertz QPO have been detected in 
three sources (Jonker, Mendez, & van der Klis 2000). 

Kilohertz QPOs have been detected with high confidence in five sources 
that also produce burst oscillations (see Table 1). In all five, the burst oscilla- 
tion frequency v osc is within 4% of the maximum observed value of Av or within 
0.35%-15% of 2Av 

max- There are no sources in which Av has been measured 
with high confidence and found to be incommensurable with v osc . Commensu- 
rability of v osc and Av appears to be a general feature of these oscillations. 

Low-frequency QPOs. — Three types of QPOs or power spectral bumps with 
frequencies <j 100 Hz have been detected in many neutron star sources (see van 
der Klis 1989, 2000): (1) the ~6-20 Hz normal- and flaring-branch oscillations 
observed in power spectra of the X-ray emission from the so-called Z sources 
when they are in their normal and flaring branch spectral states, (2) the ~ 15- 
60 Hz horizontal branch QPO seen in the Z sources when they are in their 
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Table 


1. 


Frequencies of 


Burst Oscillations 


and Kilohertz QPOs a 


Source 




^osc^max (Hz) 


A^ ma x (Hz) 


^osc/^^max 


References 


4U 1702- 


43 


330.55 ± 0.02 


344 ±7 


0.96 ± 0.02 


1 


4U 1728- 


34 


364.23 ± 0.05 


349.3 ± 1.7 


1.043 ±0.005 


2 


4U 1608- 


52 


619 


301.3 ±7.9 


2.05 ±0.013 


3,4,5 


4U 1636- 


53 


581.75 ±0.13 


254 ±5 


2.29 ± 0.04 


6 






290 b 




1.14 


7 


KS 1731- 


■26 


523.92 ± 0.05 


260 ± 10 


2.015 ±0.077 


8 



"Updated from van der Klis (2000). 6 The oscillation at 290 Hz appears to be 
intermittent (see Strohmayer 2001a, 2001b); when present, its amplitude is 
~40% of the amplitude of the oscillation at 580 Hz. References: [1] Mark- 
wardt et al. (1999) [2] Mendez & van der Klis (1999) [3] Chakrabarty et al. 
(2000) [4] Mendez et al. (1997) [5] Mendez (2000, personal communication) 
[6] Mendez et al. (1998) [7] Miller (1999) [8] Wijnands & van der Klis (1997) 



horizontal-branch spectral state, and (3) the ~ 20-80 Hz peaked noise compo- 
nents and low-frequency QPOs (LFQPOs) observed in the atoll sources. The 
latter two types of QPOs may be related to one another. The frequencies of 
these QPOs and other power spectral features are positively correlated with the 
X-ray flux and the frequencies of the kilohertz QPOs (van Straaten et al. 2001). 

2.2. Brightness Oscillations in Black Hole Candidates 

A variety of QPOs have been detected in black hole candidates (see Remillard 
et al. 2001; van der Klis 1998, 2000). QPOs with frequencies ^50 Hz have 
been detected in half a dozen black hole candidates. In addition, power spectra 
of the X-ray brightness of several black hole candidates and the peculiar source 
Cir X-l show two QPOs or broad noise peaks with frequencies between ~ 0.1 Hz 
and 10 Hz. 

High-frequency QPOs. — Pairs of high-frequency QPOs with relatively stable 
frequencies have been detected in two black-hole systems (see Table 2). QPOs 
at 40 Hz (Strohmayer 2001d) and 67 Hz (Morgan, Remillard, & Greiner 1997; 
Remillard Sz Morgan 1998) have been seen in GRS 1915±105. Their frequencies 
remain almost constant as the X-ray flux changes (any change in the frequency 
of the 40 Hz oscillation during the sequence of observations must have been 
<^ 1%; the frequency of the 67 Hz oscillation varies by only a few percent as the 
X-ray flux varies by a factor of several). Stable QPOs at 300 Hz (Remillard et 
al. 1999) and 450 Hz (Strohmayer 2001c) have been observed in GRO J1655— 40. 

Single, relatively high-frequency QPOs or peaked noise features have been 
detected in XTE J1550-564 (Remillard et al. 1999a; Homan et al. 2001; Miller 
et al. 2001), 4U 1630-47 (Remillard & Morgan 1999), XTE J1859±226 (Cui et 
al. 2000), and Cyg X-l (Remillard et al. 1999a, 1999b). The frequencies of some 
of these high-frequency QPOs are known to vary substantially; others have not 
been seen to vary but have been observed only briefly. 
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Table 2. Hig 


h-Frequency 


QPO Pairs in 


Black Hole 


Candidat 


es 


Source 


v (Hz) rms amplitude 


f/FWHM 


Av/u a 


Ref. 


GRS 1915+105 


40 


£2% 


~ 5-8 


£l% 


1 




67 


~ 1-2% 


~ 6-20 


~ 5% 


2,3 


GRO J1655-40 


300 


~ 0.8% 


~ 4 


^10% 


4 




450 


-5% 


23 


^4% 


5 



a Av is the observed variation in the QPO centroid frequency v or the up- 
per limit on this variation. References: [1] Strohmayer (2001d) [2] Mor- 
gan et al. 1997 [3] Remillard & Morgan 1998 [4] Remillard et al. 1999 
[5] Strohmayer (2001c) 



Low-frequency QPOs. — A variety low-frequency QPOs or peaked noise com- 
ponents have been identified in the power spectra of black hole candidates (see 
Psaltis, Belloni, k van der Klis 1999; van der Klis 2000). 

2.3. High-Frequency QPO Pairs in Neutron Stars and Black Holes 

The discovery of pairs of high-frequency QPOs in the two black hole sources 
GRO J1655— 40 and GRS 1915+105 raises the question whether they are pro- 
duced by the same physical mechanism as the kilohertz QPO pairs seen in more 
than twenty neutron star sources (Strohmayer 2001d). The frequencies of these 
QPOs are roughly in the range expected for dynamical frequencies in the vicinity 
of stellar-mass black holes and ~ 2Mq neutron stars, respectively. 

One indication that the high-frequency QPO pairs seen in the black holes 
have a different origin from the kilohertz QPO pairs observed in neutron stars 
is that their phenomenology is significantly different (see van der Klis 2000; 
Strohmayer 2001c, 2001d). Perhaps most significantly, the frequencies of the 
kilohertz QPOs increase by factors ~ 2-3 in some sources when the X-ray flux 
increases. In contrast, the frequencies of the high-frequency black hole QPOs 
do not change much if at all when the X-ray flux increases. Pairs of kilohertz 
QPOs are seen in both the Z and atoll sources, even though the latter sources 
are up to ~ 300 times less luminous than the former. In contrast, the 300 Hz 
QPO in GRO J1655— 40 is detected when the source is flaring. 

The rms relative amplitudes of the kilohertz QPOs in the 2-60 keV energy 
band are sometimes much greater than their FWHM/Vqpo ratios. This excludes 
models in which the QPO is produced by oscillations in the luminosity of gas in a 
ring in the disk that is narrow enough so that the range of orbital frequencies in 
the ring does not exceed the FWHM of the QPO. In contrast, the FWHM/i/qpo 
ratios of the high-frequency black hole QPOs are generally much larger than 
those of the kilohertz QPOs and all exceed the rms relative amplitude of the 
respective QPOs. Finally, the rms relative amplitudes of both kilohertz QPOs 
increase with X-ray energy up to at least ~ 12-15 keV. The rms amplitude of 
the 300 Hz QPO in GRO J1655— 40 is stronger in the soft X-ray band whereas 
the 450 Hz QPO is detected only above ~ 13 keV; conversely, the 67 Hz QPO 
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in GRS 1915+105 is detected in the soft X-ray band, whereas the 40 Hz QPO 
is only detected at energies above ~ 13 keV. The question of whether the high- 
frequency QPO pairs seen in neutron star and black holes sources are physically 
related is discussed further in §3.4. This issue clearly deserves further attention. 

3. QPO Models and Comparisons with Observation 

A variety of models have been proposed to explain the various QPOs observed 
in neutron star and black hole LMXBs. In this section I describe several QPO 
mechanisms that have been proposed and compare them with some of the data 
on QPOs and other power spectral features in neutron star and black hole 
sources. I then discuss briefly our general understanding of QPOs in neutron 
stars and black holes and their relationships to one another. 

3.1. QPO Models 

Sonic-point beat-frequency (SPBF) model. — General relativistic gas dynamical 
and radiation transport calculations (Miller & Lamb 1993, 1996; see also Tsuribe 
et al. 1995 and references therein) show that at low accretion rates, the azimuthal 
drag force exerted by radiation from a relativistic star on gas orbiting near 
it removes angular momentum from the gas sufficiently quickly that the gas 
accelerates radially inward, causing the radial flow to become supersonic before 
it reaches the star or the ISCO, if one exists. This is an effect of strong-field 
general relativity (a transition of this type is not possible in Newtonian gravity). 
To lowest order, the radius r sp of the sonic point is determined by balancing 
the angular momentum brought into the interaction region by the accretion 
flow, which is oc M£ acc , where £ acc is the mean specific angular momentum of 
the accreting gas, and the angular momentum removed from the flow by the 
radiation drag, which is cx L acc oc M. Hence to lowest order, r sp is independent 
of M. However, the optical depth of the flow from the disk to the stellar surface 
increases with M, reducing the radiation drag as M increases. This causes r sp to 
decrease with increasing M. At sufficiently high M, the accretion flow remains 
nearly circular until it reaches the ISCO, if one exists. Once it reaches the 
ISCO, gas accelerates inward and becomes supersonic even if it does not lose 
any angular momentum to the radiation field and hence r sp ~ r; sco = const. 
The existence of a sonic point is an effect of the strong gravitational field when 
f S p ~ nsco as well as when r sp > ri sco . If the neutron star is so large that there 
is no ISCO, the region of nearly circular flow may extend to the stellar surface. 

Miller, Lamb, k Psaltis (1998b, 1998c) and Lamb, Miller, & Psaltis (1998b) 
showed that inhomogeneities in the accretion flow near r sp produce two related 
high-frequency QPOs. Where gas falling inward from clumps near r sp impacts 
the stellar surface, it produces bright footprints that move around the star with 
a frequency nearly equal to the general relativistic orbital frequency f or b of the 
clumps at r sp , causing the X-ray flux seen by distant observers to oscillate as the 
footprint moves into and out of view. The frequency of this oscillation is usually 
less than the frequency ^xc(^sp) of a general relativistic Keplerian circular orbit 
at r sp , because the radial component of the radiation force causes f or b to be less 
than VKc{ r sp) and the radially inward drift of the clumps causes the footprint 
to move around the star with a frequency less than f or b (Lamb & Miller 2001). 



8 



Frederick K. Lamb 



The weak magnetic field of the star funnels extra gas toward the magnetic 
poles, producing a broad beam of slightly stronger radiation that rotates with 
the star. For disk flows that are prograde relative to the star's spin, clumps 
orbiting near r sp move through this beam with frequency vb = f or b — v s . When 
a clump is illuminated by the beam, gas in it loses angular momentum to the 
radiation at a faster rate, so more gas falls inward, causing the luminosity of the 
footprint to oscillate with a frequency close to vb- The inward radial velocity of 
the flow usually makes the frequency slightly greater than vb (Lamb & Miller 
2001). If the disk flow near the star is temporarily retrograde, clumps orbiting 
near r sp move through the beam with frequency v' B = f or b + v s , causing the 
luminosity to oscillate with a frequency close to but usually greater than v' B . 

Light travel time effects are negligible and hence the time between two 
successive maxima of the mass inflow rate from a given clump is very close to 
the beat period At-Q = l/(^orb — v s ). Let t\ and t 2 = t\ + At-Q be the times of 
successive maxima in the mass inflow rate. Then the frequencies of the lower 
and upper kilohertz QPOs are 

1 02 " 01 + (A0 2 - A&) 

v \ = n — r an d v 2 = n 1 — r - > 

t 2 - h + (At 2 - Aii) t 2 - h + {At 2 - Ati) 

where 2-7T01 and 2ir(j) 2 are the azimuthal coordinates of the clump at t\ and t 2 , 
Ati and At 2 are the times required for the gas stripped from the clump at t\ 
and t 2 to reach the surface of the star, and 2nA(j) l and 2-7rA0 2 are the changes 
in the azimuthal phase of the gas released at t\ and t 2 as it falls from the clump 
to the stellar surface. 

The inward drift of a given clump causes the time required for gas to spiral 
inward from it to the stellar surface to decrease steadily. Hence At 2 < At± . Disk 
flow is usually expected to be prograde relative to the star's spin. Then, to first 
order in the inward drift velocity v c \ of a clump, v\ is z^b(1 + v c \d r At) > vq', v\ 
is greater than v& because the inward motion of the clump reduces the spatial 
separation between successive maxima of the mass flow rate, analogous to the 
Doppler shift of a sound wave. Numerical simulations (Lamb Sz Miller 2001) 
show that v 2 is usually less than f rb 5 because the gas falling inward from the 
clump winds a progressively smaller distance around the star as the clump drifts 
inward. Hence Av is less than u s . If instead the accretion flow near the star is 
temporarily retrograde, v\ is expected to be comparable to but less than f rb; 
while v 2 is usually comparable to but greater than z/b p = v OI ^ + v s . Hence in 
this case Av is usually greater than v s . 

In the SPBF model, the frequencies of the X-ray oscillations depend on the 
mass M and angular momentum J of the star, the EOS of neutron star matter, 
the magnetic field of the star, and M. Numerical simulations (Miller et al. 1998b; 
Lamb & Miller 2001) indicate that in its simplest form, the SPBF model gives 
r sp ~ l-3-R star and X-ray oscillation frequencies ~ 500-1500 Hz. Deviations of 
v\ and v 2 from vb and f or b by <^ 5% caused by the inward motion of the accretion 
flow at r sp are sufficient to account for the observed differences between Av and 
v s and the decrease of Av with increasing v 2 described in §2.1. The X-ray flux 
for both kilohertz QPOs is produced chiefly at the stellar surface, where most 
of the gravitational energy is released, and hence both can have relatively high 
rms amplitudes. The existence of a sonic point in the flow is a strong-field 
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general relativistic effect, and the frequencies of the sonic-point kilohertz QPOs 
are governed primarily by the spacetime near the neutron star. Hence they are 
sensitive probes of this spacetime, including whether there is an innermost stable 
circular orbit (ISCO) and the gravitomagnetic torque created by the spin of the 
star. In the SPBF model, the low-frequency QPOs are assumed to be produced 
by the magnetospheric beat-frequency mechanism described below. 

Relativistic precession (RP) models. — These models identify the frequencies 
of QPOs with the basic frequencies of free-particle geodesic motion around the 
compact object. The frequencies considered are the general relativistic mean 
azimuthal (Keplerian) frequency vk, radial epicyclic frequency v r , and merid- 
ional oscillation frequency v^. For circular geodesies around a Kerr black hole 
of mass M, these frequencies (in units of c^/GM) are i>k (M,j,r) = [2vrf 3 / 2 (l + 
jr _3/2 )] _1 , O r (M,j,r) = \vK\il-Qf- 1 + 8jf~ 3 / 2 -^f' 2 ) 1 / 2 , and 0^(M,j,r) = 
\vk\(1 ~~ 4jf -3 / 2 + Sj 2 ?^ 2 ) 1 / 2 , where f is the Boyer-Lindquist radius in units 
of GM/c 2 and j = cJ/GM 2 (see Markovic 2000). The frequencies of eccen- 
tric geodesies around black holes cannot be expressed analytically. For neutron 
stars, even the frequencies of circular geodesies must be computed numerically, 
using general relativistic stellar models. 

In the clump precession model (Stella & Vietri 1998, 1999; Cui, Zhang, & 
Chen 1998b; Karas 1999; Stella, Vietri, & Morsink 1999), the frequencies of 
QPOs and other features in the power spectra of neutron stars and black holes 
are the frequencies of purely geodesic motion of gas clumps around the compact 
object. The upper and lower kilohertz QPO frequencies are assumed to be the 
Keplerian frequency uk and the apsidal precession frequency uap = V K — v r \ 
the frequency of any low-frequency QPO is assumed to be the nodal precession 
frequency z^np = vk — Vty or 2^np- These frequencies depend on M, j, and the 
periastron and apastron radii r p and r a of the geodesic; for neutron stars, they 
also depend on the EOS of neutron star matter. 

In the ring precession model (Psaltis & Norman 2000) , gas in a narrow ring 
at a special radius in the disk is assumed to behave as a separate entity. A ring 
with the properties proposed has a dense spectrum of resonances, including res- 
onances near the frequencies 2z/np> ^ap 5 and vk of an infinitesimally eccentric 
geodesic at the radius of the ring (motion along any geodesic with a finite eccen- 
tricity is quickly damped; see Markovic & Lamb 2000). In the ring precession 
model, these three frequencies are assumed to be the dominant X-ray oscillation 
frequencies and the frequencies of the low-frequency and kilohertz QPOs. These 
frequencies depend on M, j, and the radius r of the nearly circular geodesic; for 
neutron stars, they also depend on the EOS of neutron star matter. 

At present, the physics that determines which geodesies are populated by 
clumps in the precessing clump model or the radius of the ring in the precessing 
ring model has not been specified. Hence, in their present form these models 
make no definite predictions of the frequencies that should be seen or whether 
they depend on the accretion rate and luminosity. 

Disk oscillation models. — These "diskoseismic" models associate QPOs with 
the oscillation modes of the inner accretion disk (see Kato & Fukue 1980; Kato 
1990, 1993; Ipser & Lindblom 1991, 1992; Kato & Honma 1991; Ipser 1994, 
1996; for reviews, see Kato, Fukue, & Mineshige 1998; Markovic & Lamb 1998; 
Wagoner 1998). Much of this work has focused on the low-order modes of 
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geometrically thin, pseudo-Newtonian or relativistic (Novikov Sz Thorne 1973; 
Page & Thorne 1974) disks around black holes, but there has been some work on 
disks around neutron stars. The frequencies of the low-order g, c, and p-modes 
depend on the mass and angular momentum of the compact object but hardly 
at all on the accretion rate, because they are trapped waves that exist only near 
the ISCO. In studies of these modes, viscous energy dissipation is considered in 
computing the vertical structure of the disk, but angular momentum transfer 
by the viscous shear stress has usually been neglected. In this approximation, 
the motion of the fluid in the unperturbed disk is circular, i.e., there is no 
accretion. The tt;-modes of the inner disk have been studied by treating the disk 
as a succession of rings that exchange angular momentum via advection and 
the viscous shear stress. The frequencies of the tu-modes depend on where the 
nearly circular Keplerian flow ends and can vary substantially. 

The gravity-inertial modes of geometrically thin, relativistic disks are stand- 
ing waves trapped in the region of the disk where their frequencies are < v r (see 
Wagoner, Silbergleit, & Ortega- Rodriguez 2001). This trapping is a general 
relativistic effect. Hence g- modes exist only close to the compact object. The 
axisymmetric (m = 0) g-modes are centered at the radius f g where v r reaches its 
maximum value and extend over a radial interval Af g m 3.3 (L/Le) 1 ^ 2 , where 
Le is the Eddington critical luminosity; f g depends on j, decreasing from 8 for 
j = to s=s 3 for j = 0.95. The restoring force that causes oscillatory motion 
is primarily the net gravitational-centrifugal force (see Perez et al. 1997). The 
frequencies v g of the g-modes that have only a few nodes in the vertical and 
radial directions are close to v r (r g ), depend only weakly on the radial (n) and 
vertical (I) mode numbers, and decrease by ^ 5% when the luminosity doubles 
(Wagoner et al. 2001). 

The fundamental (\m\ =1 = 1) corrugation modes of geometrically thin, 
relativistic disks are trapped waves that displace the disk in the vertical direction 
and rotate around the compact object (see Silbergleit, Wagoner, Sz Ortega- 
Rodriguez 2001; Wagoner et al. 2001). Trapping of these waves is also a general 
relativistic effect. The fundamental c- modes are trapped between rj sco and the 
outer trapping radius r c . For n = 0, M ~ IOMq, and L ~ 0.1 Le, r c « 
Hsco + 0.09j _0 ' 8 (l — j) ' 2 - Hence, only if the compact object is slowly rotating 
do these modes extend over a sizeable radial interval. In the low-j limit, the 
frequency of the fundamental c-modes is equal to the gravitomagnetic frequency 
^gm(^) = vx(r) — ^v( r ) = 47T?Vif (r) 2 at the outer boundary of the trapping 
region, i.e., v c = grn (r c ); v c decreases slightly with increasing luminosity. 

Geometrically thin, viscous accretion disks around spinning compact objects 
have a family of high-frequency warping modes that are localized near the inner 
edge of the near-Keplerian flow (Markovic & Lamb 1998; see also Armitage & 
Natarajan 1999). These w- modes are spiral corrugations that rotate around the 
compact object with a frequency close to the gravitomagnetic (frame-dragging) 
frequency z> gm at the radius where their amplitudes are largest. The shapes 
of these u>-modes are determined by competition between the gravitomagnetic 
and viscous torques acting on the gas in the disk. For black holes, the high- 
frequency w-modes are situated just outside the ISCO. For neutron stars, they 
may be situated just outside the stellar surface, ISCO, sonic radius, or mag- 
netospheric boundary, depending on which effect terminates the nearly circular 
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Keplerian flow. The highest-frequency w-modes in this family are closely spaced 
in frequency and weakly damped (Q- values up to ~ 50). They are distinct from 
but related to the low-frequency u>-modes, which only perturb the accretion disk 
far from the central compact object (see Maloney, Begelman, & Pringle 1996; 
Maloney, Begelman, & Nowak 1998; Pringle 1996; Markovic & Lamb 1998). 

Accretion disks also have high-frequency p-modes, which are restored pri- 
marily by pressure-gradient forces. The fundamental (|m| = I = 0) p-modes 
span only a very small radial interval just outside rj sco and have frequencies 
\i> p \ > v r (Wagoner et al. 2001). Their frequencies increase by ~5% when the 
luminosity doubles. 

There are a variety of mechanisms by which these disk modes could modu- 
late the X-ray flux of an accreting neutron star or black hole (see, e.g., Nowak 
k Wagoner 1993; Nowak et al. 1997). 

Magneto spheric beat- frequency (MBF) model. — This model generates a rel- 
atively coherent QPO by the interaction of the stellar magnetic field with gas in 
the accretion disk (Alpar & Shaham 1985; Lamb et al. 1985; Lamb 1991; Miller 
et al. 1998b). The frequency of the QPO is the rate at which the orientation 
of the stellar magnetic field repeats, relative to the inhomogeneities in the disk 
at the coupling radius. This frequency is the magnetospheric beat frequency 
— u s ), where f rb,m is the orbital frequency at the radius r m 
at which the gas in the disk couples strongly to the stellar magnetic field and 
k = 1, 2, . . . describes the symmetry of the magnetic field there (e.g., for an offset 
dipole, both the k = 1 and k = 2 beat frequencies would be generated). The 
frequency of the QPO depends on M, j, the neutron star matter equation of 
state (EOS), M, the disk structure, and the important multipole moments of 
the stellar magnetic field at r m . 

Coronal flow oscillation model. — Spectral modeling suggests that a central 
region of radially inflowing gas forms around neutron stars in LMXBs when their 
accretion rates approach the Eddington critical rate Me (Lamb 1989, 1991). 
Numerical studies of such flows show that they become overstable when the 
accretion rate reaches a critical value somewhat less than Me (Fortner, Lamb, 
Sz Miller 1989; Fortner 1992). The resulting oscillation modulates the X-ray 
spectrum at the frequency z/ OSCiC ~ 2/tj ~ 6 Hz, where tf is the inflow time 
from the outer edge of the radially inflowing gas. This edge is determined by 
the interaction of the radiation from the neutron star with the disk corona. 



3.2. Comparison of Models with Neutron Star QPOs 

As noted in §2.1, a key feature of all the neutron star QPOs is that their frequen- 
cies vary significantly and are positively correlated with each other and with the 
inferred X-ray flux. This imposes very significant constraints on models of the 
neutron star QPOs. The frequency of the upper kilohertz QPO is comparable 
to the frequencies of orbits near the IS CO or the surface of a 1.5-2Mq neutron 
star. Most models assume that V2 is the orbital frequency f or b of gas at some 
special radius in the accretion disk (see Lamb et al. 1998b; van der Klis 1998, 
2000). A fundamental question in such models is, what physics determines this 
special radius? Whatever it is, it must cause the radius to vary by a factor ~ 1.5 
as the luminosity of the source varies by ~ 50%. 
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Sonic-point beat- frequency model. — In the SPBF model, the special radius 
in the disk is the radius r sp of the sonic transition produced either by interaction 
of the radiation from the neutron star with the accretion flow or by the approach 
of the flow to the ISCO (see §3.1). The model predicts that the upper and lower 
kilohertz QPOs usually have frequencies close to forbtjsp) an d forbfjsp) — v s -> 
respectively. The strongest motivation for the SPBF model is the evidence that 
the burst oscillations have frequencies close to the stellar spin frequency or its 
second harmonic (see §2.1) and that the frequency separation Av = v<i — v\ 
is close to v osc or 0.5f osc (see Table 1). The steady accumulation of sources 
in which Av and v osc are commensurable and the absence of any that are not 
makes it seem unlikely that this commensurability is merely a coincidence. 

As the accretion rate through the inner disk increases, r sp decreases, 
causing both kilohertz QPO frequencies to increase, in agreement with their 
observed behavior. In addition to explaining the observed frequency ranges 
of v\ and V2 and their increase with inferred M^, the SPBF model appears 
able to account for the relatively high amplitudes and coherence of the kilohertz 
QPOs and the steep increase of QPO amplitude with photon energy (Miller et al. 
1998b). The amplitudes of the kilohertz QPOs and the amplitude and frequency 
behavior of the HBOs suggest that the latter are produced by the magnetospheric 
beat-frequency mechanism (see Miller et al. 1998b; Psaltis et al. 1999b). The 
neutron star magnetic fields indicated are comparable to those expected if the 
kilohertz QPO sources are the progenitors of the recycled millisecond pulsars 
(Miller et al. 1998b; Lorimer 2000). 

Stella & Vietri (1998) suggested that the kilohertz QPOs are generated by 
the SPBF mechanism but that the HBOs and low-frequency QPOs are produced 
by nodal precession of gas at r sp . However, further study has shown that the 
frequencies and frequency behaviors of the HBOs and LFQPOs are inconsistent 
with nodal precession, unless the spin rates of these stars are ~ 800-900 Hz, 
much larger than the rates indicated by the burst oscillations and inconsistent 
with the SPBF model (Markovic & Lamb 1998; Psaltis et al. 1999; Morsink & 
Stella 1999; Kalogera & Psaltis 1999). 

As explained in §2.1, the SPBF model predicts that Av is usually somewhat 
less than v s and decreases with increasing (in Miller et al. 1998b the inward 
radial velocity of the disk flow was neglected, in which case Av = v s ). The 
inflow velocities expected at r sp naturally produce shifts in v\ and v<i of the size 
(0.2%-5%) needed to explain the differences observed between Av and v osc an d 
the variation of Av with v<i- Figure 1 shows that fits of the SPBF model agree 
fairly well with the observed behavior of the kilohertz QPOs in the four sources 
shown. 

One puzzling aspect of the kilohertz QPOs is that v\ and V2 vary differ- 
ently with countrate (or X-ray flux) on different timescales. On short timescales 
(typically ~ hours), the frequencies of the kilohertz QPOs are tightly correlated 
with the apparent X-ray luminosity, producing relatively narrow, sloping tracks 
in plots of ^qpo vs. flux, but on longer timescales (typically ~ days), these tracks 
drift in flux by as much as 50% (see the left panel of Fig. 2). Understanding 
the cause of these drifts is important in its own right, but it is essential in or- 
der to determine whether the plateaus seen in these tracks are evidence for the 
existence of an ISCO (see §4). 
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Figure 1. Fits of the sonic-point beat-frequency gas dynamical model 
to the kilohertz QPO frequencies observed in four sources, assuming 
the disk flow is prograde relative to the neutron star; the spacetime 
was approximated by the Schwarzschild spacetime. The best-fitting 
neutron star masses and spin rates and the resulting \ 2 values are 
listed in Table 3 below. After of Lamb & Miller (2001). 



In the sonic-point model, r sp is determined by a balance between the in- 
ward flux of angular momentum through the inner disk, which is proportional 
to Mrf, and the radiation drag on the gas there, which is proportional to the lu- 
minosity. Using a phenomenological model, van der Klis (2001) has shown that 
this mechanism may be able to explain the drift of the v\ and V2 vs. countrate 
tracks, if the luminosity — and hence the radiation drag — has a component that 
is proportional to the average of the future mass flux through the inner disk, 
with an averaging timescale of about a day, as might occur if part of the mass 
flux onto the star comes from rapid infall of gas from the disk at radii > r sp (see 
the right panel of Fig. 2). This possibility clearly deserves further study. 

Probably the most significant issue for the SPBF model is the lack of com- 
pelling evidence for dynamically significant magnetic fields in these neutron 
stars. There is a wealth of indirect evidence for such fields: disk accretion 
would spin stars up to spin rates much greater than ~ 350 Hz if their magnetic 
fields were negligible (Ghosh & Lamb 1992); it is difficult to understand the 
high stability of the oscillations in burst tails if there is no marker on the sur- 
face of the star; the requirement that many neutron stars in LMXBs must have 
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Figure 2. Left: v\ vs. countrate tracks observed in 4U 11608—52 
(Mendez et al. 1999). Right: Simulated tracks (van der Klis 2001). 



dipole magnetic fields ~ 10 -10 G, if they are the progenitors of the recycled 
millisecond pulsars (Kulkarni &Phinney 1994; Lorimer 2001); the correlation of 
X-ray spectra and HBO, LFQPO, and kilohertz QPO properties with inferred 
magnetic field strength (see Psaltis, Lamb, & Miller 1995, 1998a; Miller et al. 
1998b; Lamb et al. 1998b, Fig. 2); and other evidence. However, this evidence 
is not considered compelling by some. 

The discovery of burst oscillations in SAX J1808— 369 with a frequency 
close to its 400 Hz frequency spin frequency (see §2.1) is consistent with the 
hypothesis that X-ray burst oscillations are generated by the spin of the star. 
However, no kilohertz QPOs were detected during the brief interval in 1998 
when SAX J1808— 369 was observed by RXTE, so its relation to the kilohertz 
QPO sources remains unclear at present. Continuous observations of this source 
should be a high priority when it next appears. Another issue for the SPBF 
model is the apparent similarity of the correlations between ^lfqpo and v\ in 
neutron stars and between the frequencies of some of the low-frequency band- 
limited noise components seen in the power spectra of black holes (see §3.4). 

Magneto spheric and magneto spheric beat-frequency models. — The possibil- 
ity that the frequencies of the kilohertz QPOs might be the Keplerian frequency 
z^Km at the radius r m where the gas in the disk strongly couples to the star's 
magnetic field and the magnetospheric beat frequency vg m = ^ m — v s has 
been suggested by Strohmayer et al. (1996) and Ford et al. (1997). Cui et al. 
(1998), Cui (2000), and Campana (2000) have also discussed the possibility that 
QPOs are generated at r m . However, it is difficult to understand the production 
of two fairly coherent oscillations, the large observed variation of the kilohertz 
QPO frequencies, and other properties of the kilohertz QPOs in this model (see 
Miller et al. 1998b; Lamb et al. 1998a, 1998b). 

Relativistic precession models. — In the clump precession model, the kilo- 
hertz QPO frequencies V2 and v\ are assumed to be the azimuthal and apsidal 
precession frequencies vy. and z^ap of gas clumps on slightly eccentric, tilted 
geodesies, while the frequencies of the LFQPO and HBO are the first or sec- 
ond harmonics of the nodal precession frequencies z^np of these geodesies. As 
noted in §3.1, the physics that determines how the clumps are formed and which 
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Figure 3. Upper left: Kilohertz QPO frequencies measured in 
Sco X-l (data) compared with the best-fitting vap vs. uk relations 
given by innnitesimally eccentric geodesies around a nonrotating neu- 
tron star (solid curve), a star with v s = 450 Hz constructed using the 
UU equation of state (dashed curve), and a black hole (dotted curve). 
Lower left: Circumferential radii of the geodesies used. Upper right: 
Comparison of the ul vs. v\ behavior seen in 4U 1728—34 with the 
2z^np vs. uk relations given by the geodesies that best fit the v\ vs. 
V2 data. Lower right: Apastron radii of the highly eccentric geodesies 
required to fit the v\ vs. vi data. From Markovic & Lamb (2000). 



geodesies are populated has not been specified. Hence the model makes no def- 
inite predictions about whether the frequencies of the QPOs depend on the 
accretion rate. 

The frequencies z^k and z^ap for orbits near neutron stars are in same range 
as the frequencies of the kilohertz QPOs, and it was initially thought that the 
frequencies of the kilohertz QPOs could be fit by vap^k relations constructed 
using slightly eccentric geodesies (Stella & Vietri 1999; see also Karas 1999). 
However, the ^ap^^k relations initially used were inaccurate (see Markovic & 
Lamb 2000). When v^p and uk are computed correctly, the best-fitting fre- 
quency relations that can be constructed using innnitesimally or even moder- 
ately eccentric geodesies miss many of the measured frequency pairs by ~ 30% 
(see the left panel of Fig. 3) . Only by using highly eccentric geodesies with apas- 
tron radii ~ 3-4 times their periastron radii can the kilohertz QPO frequencies 
in Sco X-l and the other kilohertz QPO sources be fit (Markovic & Lamb 2000). 
But the frequencies z^np and 2zv NP of these geodesies do not track the frequencies 
of the low- frequency QPOs (see the right panel of Fig. 3). 

Another issue is that numerical computations have shown that the dominant 
frequencies that would be generated by orbiting luminous or obscuring clumps 
are harmonics and sidebands of v r , not z^k, ^ap> and z/np or 2z^jp (Markovic 
& Lamb 2000). Also, any clumps are tidally disrupted in less than 100 orbits 
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unless they are tiny and have a mean density ^ 10 5 the interclump density. This 
limits the fractional modulation that could be produced by clumps to <C 1CP 6 . 
In the precessing clump model, the closeness of burst oscillation frequencies to 
Ais or 2Au is purely coincidental. 

The precessing ring model (Psaltis & Norman 2000) assumes there is a 
"resonant ring" of material at some special radius in the inner accretion disk and 
that this ring produces X-ray oscillations mainly at uk, vap, and 2i/np- In the 
inner disk, eccentric motions are quickly circularized. Hence only nearly circular 
motions can persist. The deviation of the best-fitting i/ap vs. z^k relations for 
approximately circular orbits from the measured kilohertz QPO frequencies by 
~ 50-100 Hz or more is therefore an issue for this model. It has been argued 
(Psaltis 2000) that hydrodynamic effects may shift v r and uk by the ~ 30% 
needed to bring the frequencies of nearly circular orbits in a ring of radial width 
~ 0.01 r into agreement with the data, if the height of the ring is ~ 0.2 r and 
varies with radius. However, a height this large appears inconsistent with the 
narrow ring approximation used in the precessing ring model (Markovic & Lamb 
2000). As noted in §3.1, the physics that creates the ring and determines its 
radius has not been specified, so the model makes no definite predictions about 
the frequencies that are generated or whether they vary with the accretion rate. 
In this model, the closeness of burst oscillation frequencies to Au or 2Av is again 
purely coincidental. 

Disk oscillation models. — The frequencies and frequency behaviors of the g, 
p, and c-modes of relativistic accretion disks are incompatible with the observed 
frequencies and frequency behaviors of the kilohertz QPOs. Unlike black holes, 
the dimensionless angular momentum parameter j (see §3.1) of realistic neutron 
stars is <J 0.2. Hence, if the Keplerian disk flow extends inward to the ISCO, 
the highest-frequency g and p- modes around 1.5-2 Mq neutron stars have fre- 
quencies ^ 600 Hz, too low to be consistent with the frequencies of the kilohertz 
QPOs. Their frequencies are also insensitive to the mass accretion rate (see 
§3.1), unlike the frequencies of all the QPOs seen in neutron stars. If radiation 
drag terminates the nearly circular Keplerian flow well outside the ISCO, these 
modes do not exist. 

The frequencies of the fundamental w-modes are highest when the nearly 
circular Keplerian flow extends inward to the ISCO (Markovic &: Lamb 1998). 
These highest frequencies are consistent with the highest observed frequencies 
of the LFQPOs (~ 80 Hz) only if these neutron stars have spin rates ~ 900 Hz. 
If the nearly circular Keplerian flow ends well outside the ISCO, the frequencies 
of the tu-modes for a given stellar spin rate will be correspondingly lower. 

Although the properties of these disk modes appear inconsistent with the 
properties of the neutron star QPOs detected to date, it is worth searching for 
them in power spectra, because detecting and identifying them would provide 
valuable information about the structure of the inner disk and the stellar spin 
rate. 

3.3. Comparison of Models with Black Hole QPOs 

The frequencies of the high-frequency QPO pairs seen in GRS 1915+105 and 
GRO J1655— 40 (see Table 2) appear to be relatively constant. The frequency 
of the single high-frequency QPO seen in XTE J1550— 564 varies from 185- 
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285 Hz (see §2.2). The frequencies of the lower-frequency QPOs seen in these 
and other sources vary and are positively correlated with one another (see §2.2 
and §3.4). These phenomenological differences suggest that different mechanisms 
are responsible for the different QPOs seen in the black hole sources. 

Sonic-point and beat-frequency mechanisms. — In current models of disk ac- 
cretion onto black holes, the radius of the transition from subsonic to supersonic 
radial inflow is close to the radius of the IS CO and typically varies by less than 
1%. However, these models neglect the drag force exerted by the radiation falling 
on the disk. Although the radiation pattern differs from the pattern in neutron 
star sources, dimensional estimates indicate that this drag force could affect the 
flow and might be strong enough to change the location of the sonic point. This 
possibility deserves further investigation 

The magnetic fields of black holes are weak and axisymmetric and accreting 
gas enters the event horizon without converting its kinetic energy to radiation, 
unlike accreting gas that collides with the surface of a neutron star. Hence 
the particular beat-frequency mechanism described in §3.1 cannot operate. It 
is conceivable that radiation from orbiting hot spots at one radius in the inner 
disk could modulate inflow from a larger radius, producing a modulation of the 
accretion with a frequency equal to the difference between the orbital frequencies 
at the two radii, but no specific models of this type have been developed. 

Orbital and relativistic precession models. — An appealing idea is that the 
frequencies of the highest-frequency black hole QPOs are the orbital frequencies 
of gas orbiting near the ISCO. This hypothesis is compatible only with QPOs 
that have nearly constant frequencies. This interpretation is possible for the 
450 Hz QPO in GRO J1655— 40 if its mass is equal to the mass estimated from 
spectral observations of the companion star (6.3 ± 0.5Mq; see Shahbaz et al. 
1999; Greene, Bailyn, & Orosz 2001) and it is rotating (Strohmayer 2001a); 
j £3 0.3 is required for this mass. On the other hand, this interpretation is 
possible for the 67 Hz QPO in GRS 1915+105 only if its mass is £ 33M (the 
lowest inferred mass is obtained if the black hole is nonrotating; if it is rotating, 
the required mass is higher) . This is about twice the measured mass of the black 
hole in this system, which is 14 ± 4Mq (Greiner 2001), indicating that if 67 Hz 
is an orbital frequency, the orbit is well outside the ISCO. 

The dominant frequencies in the precessing clump and ring models are as- 
sumed to be the orbital, apsidal precession, and nodal precession frequencies vk, 
i/ap, and f N p (see §3.1). The frequency of the 300 Hz QPO in GRO J1655-40 
could be the gravitomagnetic precession frequency at the ISCO (Cui et al. 1998), 
but only if j « 1; in this case the frequency of the 450 Hz QPO could not be the 
Keplerian frequency at the ISCO, which is ;> 1 kHz (Strohmayer 2001c). The 
frequencies of the 0.9 Hz, 40 Hz, and 67 Hz QPOs and peaked noise components 
seen in GRS 1915+105 are very roughly consistent with z/np, i/ap, and vk at 
the ISCO of a ~ 3OM black hole with j « 0.1 (Strohmayer 2001d; see the left 
panel of Fig. 4), but as noted above, the measured mass of this black hole is 
14±4M . 

The observed correlation between the frequencies of the low- frequency QPOs 
in some black holes is roughly consistent (within a factor of 2) with the relation 
between i/ap and vk for nearly circular geodesies around spinning black holes 
(Stella et al. 1999). 
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Figure 4. Left: The v r vs. vk (top) and z^np vs. vk (bottom) fre- 
quency relations generated by sequences of nearly circular geodesies 
around 24Mq and 32Mq black holes with j = 0.12. The diamonds 
show the proposed identifications of the difference between the 67 Hz 
and 40 Hz QPO frequencies in GRS 1915+105 with v r and the 0.9 Hz 
noise peak with i^np- From Strohmayer (2001d). Right: Dependence 
of the fundamental frequencies of three inner-disk oscillation modes on 
the spin of the black hole and, in the case of the p-modes, the torque 
on the inner edge of the disk. Here a is the dimensionless angular mo- 
mentum, which is denoted by j in this review. The two curves for each 
case are for L/Le = 0.01 and 1.0. From Wagoner et al. (2001). 



Disk oscillation models. — The frequencies of the p, g, and c-modes of the 
inner disk are nearly independent of the accretion rate and hence may be related 
to the stable, high-frequency QPOs seen in GRO J1655-40 and GRS 1915+105 
(see the right panel of Fig. 4 and Nowak & Wagoner 1991, 1993; Perez et al. 1997; 
Nowak & Lehr 1998; Wagoner 1999). The p-mode is expected to extend over 
only a few percent of the inner disk and the X-ray luminosity from this region is 
an even smaller fraction of the disk luminosity, so attention has focused on the 
g- and c-modes (see Wagoner et al. 2001). 

The 300 Hz and 450 Hz frequencies of the high-frequency QPOs in 
GRO J1655-40 agree with the frequencies of the g and c modes for M = 5.9 ± 
l.OM and j = 0.917 ± 0.024 (Wagoner et al. 2001), consistent with the 6.3 ± 
0.5Mq mass estimate from spectral observations of the companion star (Shahbaz 
et al. 1999; Greene, Bailyn, & Orosz 2001). The frequency of the fundamental 
w-mode is essentially the same as that of the fundamental c-mode (see §3.1), so 
the 450 Hz oscillation might also be generated by this mode. The high inclina- 
tion of GRO J 1655-40 makes it more likely that these nonaxisymmetric modes 
can produce X-ray oscillations with significant amplitudes. 

Identification of the 67 Hz high-frequency QPO in GRS 1915+105 with the 
g-mode gives M = 18.2 ± 3.1, consistent with the 14 ± 4Mq estimated mass of 
this black hole (Greiner 2001), if j = 0.701 ± 0.043 (Wagoner et al. 2001). For 
such a black hole, the frequency of the c-mode is about 50 Hz, somewhat higher 
than the 40 Hz frequency of the lower high-frequency QPO in GRS 1915+105. 
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3.4. Discussion 

A large number of QPOs have now been seen in accreting neutron stars and 
black holes and a variety of mechanisms have been suggested to explain them. 
An interesting question is whether the mechanisms that produce the QPOs in 
black holes and neutron stars are similar or different. 

QPO mechanisms in neutron stars and black holes. — Different QPO mech- 
anisms are possible in neutron stars and black holes, because of the physical 
differences between them. Moreover, even the same mechanism may operate 
differently because of these physical differences. 

The existence of a stellar surface and of magnetic fields tied to the sur- 
face are probably the physically most important properties of neutron stars not 
shared by black holes. These properties of neutron stars are likely to affect 
the frequencies, waveforms, and other properties of QPOs that originate within 
a few stellar radii, causing differences from those originating near black holes, 
even if the mechanism is basically the same. Depending on their masses and 
the equation of state of neutron star matter, some neutron stars may have radii 
large enough that there is no ISCO. In contrast, there is always an ISCO around 
a black hole. The spectra and angular patterns of the radiation coming from 
the surfaces of neutron stars and the inner disks around them are expected to 
be significantly different from those of the radiation coming from the inner disks 
around black holes. As discussed in §3.1, the drag force on accreting gas pro- 
duced by the radiation coming from the surface of a neutron star is dynamically 
important. The neutron stars in LMXBs are thought to be the progenitors of 
the recycled radio pulsars (see Ghosh & Lamb 1992; Kulkarni & Phinney 1994; 
Lorimer 2001) and to have ~ 10 8 -10 9 G dipole magnetic fields. There is also 
a fairly large body of indirect evidence for such fields (see §3.2). Surface fields 
^ 10 6 G have dynamical effects on the accretion flow near the stellar surface 
(Ghosh & Lamb 1992) and would be expected to affect the properties of QPOs 
that originate within a few stellar radii. 

The exterior spacetimes of neutron stars and black holes are also different 
and hence the frequencies of dynamical motions near them are different (see, 
e.g., Miller et al. 1998a). The two most important differences are their masses 
and angular momenta. For fixed j, the dynamical frequencies near neutron stars 
and black holes scale as 1/M. The masses of black holes in LMXBs are typically 
~ 5-10 times larger than the masses of neutron stars, so dynamical frequencies 
near them should be ~ 5-10 times lower for the same j. However, j is limited to 
^ 0.2 for realistic models of neutron stars but can approach unity for black holes. 
Gravitomagnetic frequencies, and hence the frequencies of the c- and lu-modes, 
can therefore be much higher near black holes. 

Because of these important differences between neutron stars and black 
holes, it would be surprising if the highest-frequency QPOs produced by these 
objects, which presumably are generated close to their surfaces and event hori- 
zons, respectively, have identical properties. 

Frequency correlations. — Psaltis et al. (1999a; see also van der Klis 2000) 
have called attention to correlations between the frequencies of some of the low- 
and medium-frequency QPOs and noise components observed in several black 
hole sources and the similarity of the slopes of these correlations to the slopes of 
some of the frequency correlations observed in neutron star sources (see Fig. 5). 
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Figure 5. Frequencies of some of the high-, medium-, and low- 
frequency QPOs and peaked noise components seen in the power 
spectra of the LMXBs listed, plotted against the frequencies of the 
medium-frequency QPOs seen in these same sources. The dotted line 
connects the high-frequency QPO pairs seen in GRO J1655— 40 and 
GRS 1915+105. Left: Frequencies in the range 0.1 Hz to 10000 Hz. 
Right: Expanded view of the 40-1000 Hz range. After Psaltis et al. 
(1999a). 



These frequency variations are thought to reflect changes in the accretion rate. 
Psaltis et al. have suggested that some of these power spectral features may have 
similar origins in neutron stars and black holes. 

Their high and relatively stable frequencies and other properties seem to 
set the high-frequency QPO pairs seen in GRO J1655-40 and GRS 1915+105 
apart from the other black hole QPOs (see Fig. 5 and §2.2). As discussed §3.3, 
the highest frequencies may be orbital frequencies or the frequencies of disk 
oscillations near the ISCO. It appears that the differences between the highest 
frequencies seen in these two sources are not due solely to their different masses. 
If so, the mechanisms involved must be at least somewhat different. 

The frequencies of the upper kilohertz QPOs also stand apart in Figure 5. 
These are thought to be the orbital frequency at the inner edge of the nearly 
circular flow around the neutron star (see §3.2). If the QPO mechanisms in 
neutron stars and black holes are indeed the same, it is somewhat puzzling that 
this fundamental frequency has not been detected in most black holes. 

The frequency correlations of the lower kilohertz and low- frequency QPOs in 
atoll sources such as 4U 1728—34 and 4U 1608—52 are separated from the main 
correlation line established by the frequency correlations of the lower kilohertz 
and horizontal-branch QPOs in the Z sources and the peculiar source Cir X-l 
by a factor ~ 2 and seem to have a much steeper slope. The frequencies of the 
medium- and low-frequency QPOs in XTE J1550-564 and GRO J1655-40 are 
separated from the main correlation line by a factor ~ 2 and do not establish a 
clear correlation. 
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These examples indicate that not all QPO frequencies are correlated in the 
same way and hence that it is possible that different mechanisms are involved in 
different types of sources. The main correlation line that is evident in Figure 5 
is established by correlations of the lower kilohertz and low-frequency QPOs in 
the Z sources and especially in the peculiar source Cir X-l and the correlations 
of the two low-frequency QPOs in several black holes, which fall approximately 
on the same line in the frequency-frequency plane. It is also worth noting that it 
is the frequencies of the two most prominent band-limited noise components in 
the power spectrum of Cyg X-l, not the QPO frequencies, that appear to follow 
roughly the main correlation seen in Figure 5 (see Nowak 2001). 

Physical origins of frequency correlations. — There are several possible ways 
that the frequency correlations evident in Figure 5 might be explained. One pos- 
sibility is that two QPOs are generated at a single radius and scale with radius 
in approximately the same way, so that when the radius changes with accretion 
rate, they trace out a one-parameter curve in the frequency-frequency plane. For 
example, at large distances from neutron stars and black holes (r ;> 10 GM/c 2 ), 

the nodal and apsidal precession frequencies satisfy z^np ~ i 71 " 1//5 (2/3) 6 / 5 (z/pA) 6/ ' 5 
(see §3.1 and Markovic 2000) and are therefore almost proportional to one an- 
other. If logz/NP is plotted versus logz^AP, the points will fall on a line with a 
slope of 1.2, close to the slope of the main correlation in Figure 5. Moreover, both 
z^np and z^pa are steep functions of radius (in relativistic units, 2-7tz>pa ~ 3 r~ 5 / 2 
while 27TZ>np ~ 2jf~ 3 ), so only a modest variation in the radius at which the 
QPOs are produced can produce a large variation in their frequencies. Such a 
mechanism might be relevant to the frequencies of the low frequency and lower- 
kilohertz QPOs generated by Cir X-l, which vary by a factor ~20. If either of 
these are orbital frequencies, the radius at which they are generated would have 
to vary by a factor ~ 7 or more, which would be difficult to understand (see be- 
low), whereas if they are precession frequencies, a variation in radius of a factor 
of 1.8 would be sufficient. The Cir X-l QPO frequencies are in approximately 
the right range to be precession frequencies (see Stella et al. 1999) and a slope 
of 1.2 may describe their correlation better than a slope of 1. 

Sometimes an apparent power-law frequency correlation over a limited range 
of frequencies may be explained by a model that does not predict a power-law 
relation. For example, the phenomenological power-law relation z^hbo v s v\ % 
accurately describes the observed correlation between the frequencies of the HBO 
and upper kilohertz QPO in the Z sources, but so does the relation ^hbo = 
v s (Kv2 X ^ — 1), which is predicted by the magnetospheric beat-frequency model 
(see Psaltis et al. 1998b). 

These considerations underscore the importance of paying attention to the 
whole spectrum of information we now have on the phenomenology of the dif- 
ferent QPOs, not just their frequencies. Of particular importance is how the 
individual QPO frequencies vary with accretion rate, because this variation is 
an important guide to the physical mechanism(s) involved. 

The challenge of variable frequencies. — I have already noted that large vari- 
ations in the frequencies of QPOs are challenging for physical theories. This is 
particularly true if the oscillations are highly coherent; understanding large vari- 
ations in the frequencies of peaked or band-limited noise components is generally 
less demanding. Several QPO models are based on the existence of one or more 
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special radii in the disk that decrease substantially as the accretion rate in- 
creases. A key question is, what physics picks out a special radius and causes it 
to vary by a substantial factor as the accretion rate changes by, e.g., ~ 50%? 

To my knowledge, the only mechanism so far identified that picks out a 
special radius that decreases markedly as the accretion rate increases is the 
radiation drag mechanism that is at the heart of the sonic-point model (see 
§3.1). In its original form, this mechanism can explain a variation in the radius 
of the sonic point by a factor ~3, from ~ R s tai to ~ 3i? s t a r- The radius and 
radius range could be greater if radiation drag only needs to act on a small 
fraction of the accretion flow, such as gas on the surface of the disk, to make a 
QPO mechanism work. It is possible that the radiation-drag mechanism could 
also work in black hole systems, although this possibility has not been studied 
in any detail. 

The radiation drag mechanism has been invoked to explain the existence 
of the special radii needed in models other than the sonic-point model. It is 
important to bear in mind that the radiation force has important dynamical 
effects on the accretion flow over an extended region, where it changes both the 
orbital and the radial velocities of the flow. These effects, as well as the existence 
of a sonic radius, must be taken into account in any self-consistent model of the 
formation of QPOs near neutron stars. 



4. Using QPOs to Study Dense Matter and Strong-Field Gravity 

As mentioned in the Introduction, almost all QPO models developed since the 
launch of RXTE have considered the effects of strong-field gravity. In this 
sense, the QPOs themselves provide information on gas dynamics and radiation 
transport in strongly curved spacetimes. In this section I discuss briefly how the 
kilohertz QPOs can be used to extract further information about the properties 
of strong gravitational fields and dense matter. 

Assuming that the upper kilohertz QPO at 1/2 is produced by orbital motion 
of gas near the neutron star, its behavior can be used to investigate orbits in 
the region of strongly curved spacetime near the star. For example, it may be 
possible to establish the existence of an innermost stable circular orbit (ISCO) 
around some neutron stars in LMXBs (see Kaaret & Ford 1997; Miller et al. 
1998b, 1998c; Lamb et al. 1998b). This would be an important step forward 
in our understanding of strong- field gravity and the properties of dense matter, 
because it would be the first confirmation of a prediction of general relativity in 
the strong-field regime. 

Possible signatures of the ISCO have been discussed by Miller et al. (1998a) 
and Lamb et al. (1998a, 1998c). The sonic-point model predicts that the fre- 
quencies of both kilohertz QPOs will increase with the accretion luminosity until 
the sonic radius — which moves inward as the mass flux through the inner disk 
increases — reaches the ISCO, at which point it will become independent of the 
accretion luminosity (see the left panel of Fig. 6) . The plateau frequency is the 
orbital frequency of the ISCO. Zhang et al. (1998), Kaaret et al. (1999), and 
Bloser et al. (1999) have reported seeing behavior similar to this in 4U 1820—30 
(see Fig. 7). 
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Figure 6. Left: Example of the dependence of the upper and lower 
kilohertz QPO frequencies on the accretion luminosity found in calcula- 
tions using the sonic-point beat-frequency model. Once the sonic point 
reaches the ISCO, the curve flattens and hence the plateau frequency is 
the orbital frequency at the ISCO. Right: Mass-radius plane, showing 
the region allowed for a nonrotating neutron star with = 1330 Hz 
(see text). The star's representative point must lie in the unhatched 
region enclosed by the solid line. After Miller et al. (1998a). 



Although the observed variations of v\ and V2 with countrate (and X-ray 
flux) are very similar to the relation predicted by the sonic-point model, several 
questions must be answered before these results can be accepted as compelling 
evidence of an ISCO. For example, how robust is the predicted relation between 
QPO frequency and Ml What is the expected behavior of the kilohertz QPO 
amplitudes as the sonic point reaches the radius of the ISCO? Naively, one would 
expect the amplitude of the QPO at the beat frequency to decrease. Could the 
apparent plateau be an artifact of shifts in the v\ and vi tracks like those seen 
in other sources (see Fig. 2)1 

Using the kilohertz QPOs to determine the masses, radii, and spin rates of 
neutron stars. — Assuming only that the frequency of a given kilohertz QPO is 
an orbital frequency, one can extract constraints on the mass and radius of the 
neutron star (Miller et al. 19998; Lamb et al. 1998). These constraints can be 
determined for any stellar spin frequency (see Miller et al. 1998b; Miller, Lamb, 
& Cook 1998), but it is easier to see how they are obtained for a nonrotating 
star. I therefore consider this case here. 

The right panel of Figure 6 shows how to construct constraints on the mass 
and radius of a nonrotating neutron star, given v%, the highest orbital frequency 
observed in the source. i? or b must be greater than the stellar radius, so the star's 
representative point must lie to the left of the (dashed) cubic curve M°(i? or b) 
that relates the star's mass to the radius of orbits with frequency v\- The high 
coherence of the oscillations constrains i? or b to be greater than R ms , the radius 
of the innermost stable orbit, which means that the radius of the actual orbit 
must lie on the M°(fi or b) curve below its intersection with the (dotted) straight 
line M°(R ms ) that relates the star's mass to R ms . These requirements constrain 
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Figure 7. Kilohertz QPO frequency vs. countrate tracks in 
4U 1820-30 observed by Zhang et al. (1998) (left) and Kaaret et al. 
(1999) (right). 



the star's representative point to lie in the unhatched, pie-slice shaped region 
enclosed by the solid line. The allowed region shown is for v\ = 1330 Hz, the 
highest value of V2 observed in 4U 0614+09 (van Straaten 2000), which is also 
the highest value so far observed in any source. 

The left panel of Figure 8 shows how this allowed region compares with 
the mass-radius relations given by five representative equations of state (for a 
description of these EOS and references to the literature, see Miller et al. 1998a). 
If 4U 0614+09 were not rotating, EOS L and M would both be excluded. How- 
ever, 4U 0614+09 is undoubtedly rotating and probably has a spin frequency 
^ 350 Hz, in which case EOS M may be allowed. However, EOS L is excluded 
for any spin rate. The right panel of Figure 8 shows how the allowed region com- 
pares with the mass-radius relations given by the FPS EOS and the most recent 
realistic EOS, denoted APR (Akmal, Pandharipande, k, Ravenhall 1998). Both 
are allowed for nonrotating neutron stars, but the FPS EOS may be excluded if 
4U 0614+09 is rapidly rotating. 

If one can establish that a given QPO frequency is the frequency of the 
ISCO, for example by detecting a stable plateau in the QPO frequency-luminosity 
relation, the mass of the star can be determined. The right panel of Figure 6 
shows how this can be done for a nonrotating star. If an observed value of V2 
is the frequency of the ISCO, the radius of the orbit involved must lie on the 
diagonal dashed line as well as on the cubic curve that corresponds to f2, i.e., 
it is at the intersection of the two. This fixes the mass of the star for a given 
proposed EOS (Miller et al. 1998a, 1998b, 1998d; Lamb et al. 1998a, 1998c). As 
a specific example, if 4U 1820—30 has a spin rate of 290 Hz and the 1060 Hz 
plateau frequency observed in this source is the frequency of the innermost sta- 
ble orbit, general relativistic calculations show that the neutron star must have 
a mass close to 2.3 Mq for the few equations of state that give a solution. If 
confirmed, such a high mass would imply that the three- nucleon interaction is 
highly repulsive at high densities and probably that there is no quark matter at 
the center of even such a massive neutron star. These conclusions would be of 
fundamental importance for nuclear physics as well as for astrophysics. 
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Figure 8. The region of the M-R plane allowed for a nonrotating 
star in which a 1330 Hz orbital frequency has been observed, com- 
pared with the mass-radius relations given by five representative neu- 
tron star equations of state (left) and two recent realistic equations of 
state (right). 

Once they have been firmly established, specific QPO models will make it 
possible to derive further constraints on the properties of neutron stars. For 
example, the fits shown in Figure 1 of the sonic-point beat-frequency model 
to the kilohertz QPO behavior observed in Sco X-l, 4U 1608-52, 4U 1728-34, 
and 4U 1820—30 give the masses and spin rates listed in Table 3. The effect of 
the star's angular momentum on its structure and the spacetime were neglected, 
which is why the inferred mass of 4U 1820-30 is 2.0 Mq, rather than 2.3 Mq. 

5. Concluding Remarks 

The launch of the Rossi X-Ray Timing Explorer has brought a flood of new 
information about neutron stars and black holes in LMXBs, stimulating imag- 



Table 3. 


Inferred Properties 


of Neutron St 


ars 


Source 


Mass (M/M ) 


Spin Rate 


X 2 /dof 


Sco X-l 


1.59 


352 


95.5/46 


4U 1608-52 


1.98 


309.5 a 


20.6/10 


4U 1728-34 


1.74 


364 


8.5/6 


4U 1820-30 


2.0 6 


279 


7.0/18 



a Set equal to one-half the burst oscillation frequency. b Set 
equal to the mass given by equating the frequency of the 
plateau to the frequency of the ISCO. 
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inative new thinking about these systems. A variety of mechanisms have been 
proposed to explain the many QPOs that have been discovered. We are only 
beginning to explore and test them. It is likely that some mechanisms which 
generate QPOs have not yet been identified. New observational results that will 
help advance our understanding are being reported almost daily. 

Despite uncertainty about the precise physical mechanisms responsible for 
generating many of the QPOs seen in neutron star and black hole systems, there 
is good evidence that the upper kilohertz QPO in neutron stars is produced by 
orbital motion of gas in the strong gravitational field near the star. Making only 
this minimal assumption, important new constraints on the masses and radii 
of the neutron stars in LMXBs can be derived. There are indications that the 
signature of an innermost stable orbit has been seen in 4U 1820—30. If this is 
confirmed, it will be the first detection of a strong-field general relativistic effect. 

The ability to study variations in the emission of neutron stars and black 
holes on the dynamical time scales near them has made it possible to address 
questions of fundamental physics. It is important to continue this effort. 
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